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Abstract 

We investigate the power of the Wang tile self-assembly model at temperature 1, a threshold 
value that permits attachment between any two tiles that share even a single bond. When 
restricted to deterministic assembly in the plane, no temperature 1 assembly system has been 
shown to build a shape with a tile complexity smaller than the diameter of the shape. In 
contrast, we show that temperature 1 self-assembly in 3 dimensions, even when growth is 
restricted to at most 1 step into the third dimension, is capable of simulating a large class of 
temperature 2 systems, in turn permitting the simulation of arbitrary Turing machines and 
the assembly of n x n squares in near optimal O(logn) tile complexity. Further, we consider 
temperature 1 probabilistic assembly in 2D, and show that with a logarithmic scale up of 
tile complexity and shape scale, the same general class of temperature r = 2 systems can be 
simulated with high probability, yielding Turing machine simulation and 0(log^ n) assembly 
of n X n squares with high probability. Our results show a sharp contrast in achievable tile 
complexity at temperature 1 if either growth into the third dimension or a small probability of 
error are permitted. Motivated by applications in nanotcchnology and molecular computing, 
and the plausibility of implementing 3 dimensional self-assembly systems, our techniques may 
provide the needed power of temperature 2 systems, while at the same time avoiding the 
experimental challenges faced by those systems. 
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1 Introduction 



Self-assembly is the process by which simple objects autonomously assemble into an organized 
structure. This phenomenon is the driving force for the creation of complex biological organisms, 
and is emerging as a powerful tool for bottom up fabrication of complex nanostructures. One of 
the most fruitful classes of self-assembly systems is DNA self-assembly. The ability to synthesize 
DNA strands with specific base sequences permits a highly reliable technique for programming 
strands to assemble into specific structures. In particular, molecular building blocks or tiles can 
be assembled with distinct bonding domains I^^ ^ ^^ L These DNA tiles can be designed to simulate 
the theoretical bonding behavior of the Tile Assembly Model^. 

In the Tile Assembly Model, the particles of a self-assembly system are modeled by four sided 
Wang tiles for 2D assembly, or 6 sided Wang cubes in 3D. Each side of a tile represents a distinct 
binding domain and has a specific glue type associated with it. Starting from an initial seed tile, 
assembly takes place by attaching copies of different tile types in the system to the growing seed 
assembly one by one. The attachment is driven by the affinities of the glue types. In particular, 
a tile type may attach to a growing seed assembly if the total bonding strength from all glues 
abutting against the seed assembly exceeds some given parameter called the temperature. If the 
assembly process reaches a point when no more attachments are possible, the produced assembly 
is denoted terminal and is considered the output assembly of the system. 

Motivated by bottom up nanofabrication of complex devices and molecular computing, a num- 
ber of fundamental problems in the tile assembly model have been considered |5 | 9 | 2i | 23 | 24 | 27 | 34 ] ^ ^ 
few problems are as follow: (1) shape fabrication, given a target shape T, design a system of tile 
types that will uniquely assemble into shape T that uses as few distinct tile types as possible; 2) 
molecular computing , given some assembled input assembly that encodes a description of a 
computational problem, design a tile system that will read this input and assemble a structure 
that encodes the solution to the computational problem; (3) shape replication^^, given a single 
copy of a preassembled input shape or pattern, efficiently create a number of replicas of the input 
shape or pattern. 

While a great body of work has emerged in recent years considering problems in the tile 
assembly model, almost all of this work has focussed on temperature 2 assembly in which tiles 
require 2 separate positive strength glue bonds to attach to the growing seed assembly. This 
is in contrast to the simpler temperature 1 model which permits attachment given any positive 
strength bond. It seems that some fundamental increase in computational power and efficiency is 
achieved by making the step from temperature 1 to temperature 2. In fact, there is no known 2D 
construction to deterministically assemble a width n shape in fewer that n distinct tiles types at 
temperature 1. This is in contrast to efficient temperature 2 systems which assemble large classes 
of shapes efficiently, including n x n squares in optimal ^( io^g^g„ ) tile types. In fact, the ability 
to simulate universal computation and assemble arbitrary shapes in a number of tile types close 
to the Kolmogorov complexity of a given shape at temperature 2 in 2D ESI has resulted in limited 
interest in exploring 3D assembly, as it would appear no substantial power would be gained in the 
extra dimension. 

While temperature 2 assembly yields very efficient, powerful constructions, it is not without its 
drawbacks. One of the main hurdles preventing large scale implementation of complex self-assembly 
systems is high error rates. The primary cause of these errors in DNA self-assembly systems stems 
from the problem of insufficient attachments of tile types liI B l 22 | 33| _ xhat is, in practice tiles often 
attach with less than strength 2 bonding despite carefully specified lab settings meant to prevent 
such insufficient bonds. Inherently, this is a problem specific only to temperature 2 and higher 
systems. For this reason, development of temperature 1 theory may prove to be of great practical 
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Table 1: In this table we summarize the state of the art in achievable tile complexities and 
computational power for tile self-assembly in terms of temperature 1 versus temperature 2 assembly, 
2-dimensional versus 3-dimensional assembly, and deterministic versus probabilistic assembly. Our 
contributions are contained in rows 2 and 3. 
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interest. 

Because of a perceived lack of power, temperature 1 assembly has received little attention com- 
pared to the more powerful temperature 2 assembly. In addition, directions such as 3D assembly 
have not received substantial attention stemming from a perceived lack of ability to increase the 
functionality of the already powerful temperature 2 systems. Interestingly, we find that both di- 
rections are fruitful when considered together; temperature 1 assembly systems in 3D are nearly 
as powerful as temperature 2 systems, suggesting that both the perception of limited temperature 
1 power and the perception of limited value in 3D are not completely accurate. 

Our Results. In this paper we show that temperature 1 deterministic tile assembly systems in 
3D can simulate a large class of temperature 2 systems we call zig zag systems. We further show 
that this simulation grants both: (1) near optimal O(logn) tile type efficiency for the assembly 
of n X n squares and (2) universal computational power. Further, in the case of 2D probabilistic 
assembly, we show similar results hold by achieving O(log^n) efficient square assembly and the 
ability to efficiently simulate any time bounded Turing machine with arbitrarily small chance of 
error. The key technique utilized in our constructions is a method of limiting glue attachment by 
create geometrical barriers of growth that prevent undesired attachments from propagating. This 
technique is well known in the field of chemistry as steric hindrance or steric protection CSHlZlS] 
where a chemical reaction is slowed or stopped by the arrangement of atoms in a molecule. These 
results show that temperature 1 assembly is not as limited as it appears at first consideration, and 
perhaps such assemblies warrant more consideration in light of the potential practical advantages 
of temperature 1 self-assembly in DNA implementations. 

Practical Drawbacks of Temperature 1 Self- Assembly. While temperature 1 assembly 
avoids many of the practical hurdles limiting temperature 2 assembly, temperature 1 assembly 
also introduces new issues. In particular, the problem of multiple nucleation, in which tiles begin 
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to grow without the presence of the seed tile, is a more substantial problem at temperature 1. 
However, further research into temperature 1 assembly may suggest and motivate new techniques 
to limit such errors. In the specific case of multiple nucleation, we discuss in this paper as future 
work a new design technique for temperature 1 self-assembly to limit such errors, even in a pure 
temperature 1 assembly model. By fully exploring techniques such as this, and any new techniques 
that may emerge, temperature 1 self-assembly may emerge as a practical alternative to temperature 
2 assembly. 

Related Work Some recent work has been done in the area of proving lower bounds for tem- 
perature 1 self-assembly. Doty et. al^^ show a limit to the computational power of temperature 1 
self-assembly for pumpable systems. Munich et. al^^ show that temperature 1 assembly of a shape 
requires at least as many tile types as the diameter of the assembled shape if no mismatched glues 
are permitted. In terms of positive results, Chandran et. al^l consider the probabilistic assembly 
of lines with expected length n (at temperature 1) and achieve O(logn) tile complexity. Kao and 
SchwellerSS and DotyJ^ use a variant of probabilistic self-assembly (at temperature 2) to reduce 
distinct tile type complexity. Demaine et. al^^ and Abel et. al^ utilize steric hindrance to assist 
in the assembly and replication of shapes over a number of stages. 

Paper Layout. In Section [2] we define the Tile Assembly Model, in Section [3] we describe an 
algorithm to convert a temperature 2 zig zag system into an equivalent temperature 1 3D system, 
or a probabilistic 2D system, in Section |4] we show how temperature 1 systems can efficiently 
assemble n x n squares, in Section [5] we show that temperature 1 systems can simulate arbitrary 
Turing machines, in Section [6] we discuss preliminary experimental simulations, and in Section [7] 
we discuss further research directions. 

2 Basics 

2.1 Definitions: the Abstract Tile Assembly Model in 2 Dimensions 

To describe the tile self-assembly model, we make the following definitions. A tile type t in the 
model is a four sided Wang tile denoted by the ordered quadruple (north(t), east(t), south(t), west(t)) 
The entries of the quadruples are glue types taken from an alphabet S representing the north, east, 
south, and west edges of the Wang tile, respectively. Each pair of glue types are assigned a non- 
negative integer bonding strength (0,1, or 2 in this paper) by the glue function to {0, 1, . . .}. It 
is assumed that G(x, y) = G{y, x), and there exists a null in S such that \/x G S, G'(null, x) = 0. 
In this paper we assume the glue function is such that G{x, y) = when x ^ y and denote G{x^ x) 
by G{x). 

A tile system is an ordered triple (T, s, r) where T is a set of tiles called the tileset of the 
system, r is a positive integer called the temperature of the system and s G T is a single tile called 
the seed tile. |T| is referred to as the tile complexity of the system. In this paper we only consider 
temperature r = 1 and r = 2 systems. 

Define a configuration to be a mapping from to T |J {empty}, where empty is a special tile 
that has the null glue on each of its four edges. The shape of a configuration is defined as the 
set of positions {i,j) that do not map to the empty tile. For a configuration C, a tile i G T is 
said to be attachable at the position {i,j) if C{i,j) = empty and G(north(t), south(C(i, j + 1))) + 
G'(east(i),west(C(i + l,j))) + G(south(i), north(C(z, j - 1))) + G(west(i), east(C(i - l,j))) > r. 
For configurations G and G' such that C{x,y) = empty, G'{i,j) = G{i,j) for all (i,j) / {x,y), 
and G'{x,y) = t for some t & T, define the act of attaching tile t to C at position {x,y) as the 
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transformation from configmation C to C . For a given tile system T, if a supertile B can be 
obtained from a supertile A by the addition of a single tile we write A — t-t B. Further, we denote 
A -^T as the set of all B such that A -^j- B and — t-^ as the transitive closure of — 

Define the adjacency graph of a configuration C as follows. Let the set of vertices be the set 
of coordinates such that C{i,j) is not empty. Let there be an edge between vertices (xi,7/i) 
and (3^2, 2/2) iff l^^i — 2:21 + |yi — 2/2I = 1- We refer to a configuration whose adjacency graph is finite 
and connected as a supertile. For a supertile S, denote the number of non-empty positions (tiles) 
in the supertile by size(S'). We also note that each tile t G T can be thought of as denoting the 
unique supertile that maps position (0, 0) to t and all other positions to empty. Throughout this 
paper we will informally refer to tiles as being supertiles. 

2.2 The Assembly Process 

Deterministic Assembly Assembly takes place by growing a supertile starting with tile s at 
position (0,0). Any t gT that is attachable at some position {i,j) may attach and thus increase 
the size of the supertile. For a given tile system, any supertile that can be obtained by starting 
with the seed and attaching arbitrary attachable tiles is said to be produced. If this process comes 
to a point at which no tiles in T can be added, the resultant supertile is said to be terminally 
produced. For a given shape T, a tile system F uniquely produces shape T if for each produced 
supertile A, there exists some terminally produced supertile A' of shape T such that A — A' . 
That is, each produced supertile can be grown into a supertile of shape T. The tile complexity of 
a shape T is the minimum tile set size required to uniquely assemble T. For an assembly system 
which uniquely assembles one supertile, the system is said to be deterministic. 

Probabilistic Assembly For non-deterministic assembly systems, we define the probabilistic 
assembly model to place probability distributions on which tiles attach throughout the assembly 
process. To study this model we can think of the assembly process as a Markov chain where each 
producible supertile is a state and transitions occur with non-zero probability from supertile A to 
each B G A — )-t- For each B € A — )-t, let ts denote the tile added to A to get B. The transition 
probability from j4 to -B is defined to be 

TRANS(A,B) = JJ^^y 

The probability that a tile system T terminally assembles a supertile A is thus defined to be 
the probability that the Markov chain ends in state A. Further, the probability that a system 
terminally assembles a shape T is the probability the chain ends in a supertile state of shape T. 

2.3 Extension to 3 Dimensions 

Extending the model, we can consider assembly in 3 dimensions by considering Wang cubes with 
added "up" and "down" glues (up(t) and down(t)), and configurations of tile types as mappings 
from to a tile set T. In this model, assembly begins as before with an initial seed cube 
(informally we will refer to cubes as tiles) at the origin, with cubes attaching to the north, west, 
south, east, top, or bottom of already placed tiles if the total strength of attachment from the glue 
function meets or exceeds the temperature threshold r. 

In this paper we only consider temperature r = 1 assembly systems in 3D. Further, we only 
consider systems that assemble tiles within the z = or z = 1 plane. To depict 3 dimensional tile 
sets and 3 dimensional assemblies, we introduce some new notation in Figure [TJ 
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1 I 2 I 3 I 4 

Figure 1: For 3 dimensional tile assemblies, we create figures with the notation depicted above. 
First, large tiles denote tiles placed in the z = plane, while the smaller squares denote tiles placed 
in the z = 1 plane. The red connectors between bottom tiles denotes some unique glue shared by 
the tiles in the figure, as does the thin black line connecting tiles in the top plane. Blue circles 
denote a unique glue connecting the bottom of the top tile with the top of the tile below it. In this 
example, the tiles are numbered showing the implied order of attachment of tiles assuming the '1' 
tile is a seed tile. 

2.4 Scaled Simulations 

Our constructions in this paper relate to algorithms that generate temperature 1 tile systems that 
simulate the assembly of a given temperature 2 tile system. Formally, we define the notion of one 
tile system simulating another system. 

Definition 1. A tile system T = (T',s',t') simulates a deterministic tile system T = (T, s,r) at 
horizontal scale factor scalcx, vertical scale factor scalcy, and tile complexity factor C if: 

1. there exists an onto function SIM : R (^T' ^ T for some subset R <ZT' such that for any non 
empty tile type T(i,j) at position {i,j) in the terminal assembly ofT, there exists exactly one 
tile type r £ R such that T{k, I) = r for some k and i such that i ■ scalcx < k < {i + l)scalex 
and j ■ scalcy < ^ < (i + l)scaley, and SIM{r) = T(i,j). 

2. for each position {i,j) such that T(i,j) = empty, it is the case that for each k and i such 
that i ■ scalcx <k<{i + l)scalex and j ■ scalcy < ^ < (j + l)scaley, then T{k,£) = empty. 

3. \T'\ < C\T\ 

Informally, a system T simulates a system F if for each tile type t placed by system F, a tile 
type representing t is place by T at the same position but scaled up by the vertical and horizontal 
scale factors. In the case that the simulating system is a 3D system, all tile positions from the 3D 
system are projected on the z = plane to apply the definition. In the case that T is a probabilistic 
assembly system, then the tile placed at a given position by T is a random variable, and we are 
interested in the probability that each tile placed by F is correctly simulated by T according to 
some given assignment of representative tile types. 

3 Simulation of Temperature r = 2 Zig-Zag Systems at Tempera- 
ture r = 1 

3.1 Zig-Zag Tile Systems 

Our first result is a construction that shows a large class of temperature 2 self-assembly systems (in 
2 dimensions) can be simulated at temperature r = 1 if growth is permitted in the third dimension. 
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Figure 2: A zig-zag tile system alternates growth from left to right at each vertical level. 



Zig-zag assembly systems. An assembly system is defined to be a zig-zag system if it has the 
following two properties: 

• The assembly sequence, which specifies the order in which tile types are attached to the 
assembly along with their position of attachment, is unique. For such systems denote the 
type of the i*^ tile to be attached to the assembly as type{i), and denote the coordinate 
location of attachment of the i^^ attached tile by {x{i),y{i)). 

• For each i in the assembly sequence of a zig-zag system, if y{i — 1) is even, then either 
x{i) = x{i — 1) + 1 and y{i — 1) = or y{i) = y{i — 1) + 1. For odd y{i — 1), either 
x{i) = x{i — 1) — 1 and y{i — 1), or y{i) = y{i — 1) + 1. 

Intuitively, a zig-zag assembly is an assembly which must place tiles from left to right up to 
the point at which a first tile is place into the next row north. At this point growth must grow 
from right to left, until the next growth north once again. Thus even rows grow from left to right, 
odd rows grow from right to left, as shown in Figure [2] 

The key property zig-zag systems have is that it is not possible for a west or east strength 1 
glue to be exposed without the position directly south of the open slot already being tiled, or, if 
this is possible, the glue is redundant and can be removed to achieve the same assembly sequence. 

3.2 3D Temperature 1 Simulation of Zig-Zag Tile Systems 

In this section we show that any temperature 2 zig-zag tile system can be simulated efficiently by 
a temperature 1 3D tile system. The scale and efficiency of the simulation is dependant on the 
number of distinct strength 1 glues that occur on either the north or south edge of a tile type in 
the input tile system's tileset T. For a tileset T, let or denote the set of strength 1 glue types that 
occur on either a north or south face of some tile type in T. The following simulation is possible. 

Theorem 1. For any 2 dimensional temperature 2 zig-zag tile system T = {T, s, 2) with ax denot- 
ing the set of distinct strength 1 north/south glue types occuring in T, there exists a 3D temperature 
1 tile system that simulates T at vertical scale = 4, horizontal scale = 0(log |(Tt|), and (total) tile 
complexity 0{\T\ log|fTT|), which constitutes a 0(log|o"j'|) tile complexity scale factor. 

Proof. To simulate a temperature 2 zig-zag tileset with a temperature 1 tileset, we replace each 
tile type in the original system with a collection of tiles designed to assemble a larger macro tile 
at temperature 1. The goal will be to "fake" cooperative temperature 2 attachment of a given tile 
from the original system by utilizing geometry to force an otherwise non-deterministic growth of 
tiles to read the north input glue of a given macro tile. 

Consider an arbitrary temperature 2 zig-zag tile system F = (T, s,2) where ctt denotes the 
alphabet of strength 1 glue types that occur on either the north or south face of some tile type 
in T. Index each glue g G ax from to |cjt| — 1- For a given g G or, let b{g) denote the index 



7 



value of g written in binary. We will refer to each glue g by the new name C(,(g) when describing 
our construction. 

To prove the theorem, we construct a tileset with 0(log|fT2'|) distinct tile types for each tile 
t G T. In particular, the number of tile types will be at most |T|(12*log |(Tt| + 12) = 0{\T\ log \(Tt\)- 

To assign tile types to the new temperature 1 assembly system, take all west growing tile types 
that have an east glue of type 'x' for some strength 1 glue 'x'. If there happens to be both east 
growing and west growing tile types that have glue 'x' as an east glue, first split all such tile types 
in to two separate tiles, one for each direction. 

For the set of west growing tile types with east glue x, a collection of tile types are added to the 
simulation set as a function of x and the subset of gt glues that appear on the south face of the 
collected tile types. The tile types added are depicted in Figure [3j In the example from Figure [3] 
there are 4 distinct tile types that share an east x glue type. As these tiles are west growing tile 
types, the temperature 2 simulation places each of these tiles using the cooperative bonding of 
glue X and glue cm in the case of the right most tile type. The east and south glue types of a west 
growing tile can be thought of as input glues, which uniquely specify which tile type is placed, in 
turn specifying two output glues, glue type a to the west and glue type cm to the north in the 
case of the right most tile type. At temperature 1, we cannot directly implement this double input 
function by cooperative bonding as even a single glue type is sufficient to place a tile. Instead, we 
use glue type to encode the east input, and geometry of previously assembly tiles to encode the 
south input. 

In more detail, the tile types specified in Figure [3] (b) constitute a nondeterministic chain of 
tiles whose possible assembly paths for a binary tree of depth log of the number of distinct south 
glue types observed in the tile set being simulated. In the given example, the tree starts with 
an input glue (x, — ). This glue knows the tile to its east has a west glue of type x, but has 
no encoding of what glue type is to the south of the macro tile to be placed. This chain of tiles 
nondeterministically places either a or a 1 tile, which in turn continues growth along two separate 
possible paths, one denoted by glue type (x,0), and the other by glue type (x, 1). By explicitly 
encoding all paths of a binary prefix tree ending with leaves for each of the south glues of the input 
tile types, the decoding tiles nondeterministically pick exactly one south glue type to pair with 
the east glue type x, and output this value as a glue specifying which of the 4 tile types should be 
simulated at this position. 

Now, to eliminate the non-determinism in the decoding tiles, we ensure that the geometry 
of previous place tiles in the z = 1 plane is such that at each possible branching point in the 
binary tree chain, exactly 1 path is blocked, thus removing the non-determinism in the assembly 
as depicted in Figure [s] (d). This prebuilt geometry is guaranteed to be in place by the correct 
placement of the simulated macro tile placed south of the current macro tile. One the proper tile 
type to be simulated is decode, the 2 output values, a and cm in the case of the rightmost tile 
type of Figure [3] (a), must be propagated west and north respectively. This is accomplished by the 
collection of tile types depicted in Figure [s] (c). Now that the north and west output glues have 
been decode, this macro tile will assemble a geometry of blocking tiles to ensure that a tile using 
this north glue as a south glue input will deterministically decode the correct glue binary string. 
In particular, pairs of tiles are placed in the plane z = 2 for each bit of the output binary string. 
The pair is place to locations vertically higher for 1 bits than for bits. The next row of macro 
tiles will then be able to decode this glue type encoded in geometry by applying a binary tree of 
decoder tiles similar to those shown in Figure [s] (b) . 

The complete conversion algorithm from a temperature 2 zig-zag system to a temperature 1 
3D system has a large number of special cases. However, the example worked out in this proof 
sketch gets at the heart of the idea. The fully detailed conversion algorithm, with all cases detailed 



8 



is described in the Appendix in Section [Cj 
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Figure 3: This tile set depicts how one tile from a collection of input tile types all sharing the 
same east glue x are simulated by combining strength 1 glues bonds with geometrical blocking to 
simulate cooperative binding and correctly place the correct macro tile. 

□ 

3.3 An Example 

To further illustrate the construction for Theorem [T] we present a simple zig-zag tile set and 
corresponding 3D, temperature r = 1 simulation tile set, described in the Appendix in Figures [6] 
and El 

3.4 2D Probabilistic Temperature 1 Simulation of Zig-Zag Tile Systems 

In this section we sketch out the simulation of any temperature 2 zig-zag tile system by a 2D 
probabilistic tile system. The basic idea is similar to the 3D conversion. Each tile type for a given 
east glue x is used to generate a binary tree of nondeterministic chains of tiles to decode a southern 
input glue encoded by geometry. However, in 2D, it is not possible to block both branches of the 
binary tree due to planarity. Therefore, we only block one side or block neither side. The basic 
idea behind the decoding is depicted in Figure |4j 

The key idea is to buffer the length of the geometric blocks encoding bits by a buffing factor 
of k. In the example from Figure |4j a 2-bit string is encoded by two length 2k stretches of tiles 
(k=4 in this example), where the binary bit value 1 is represented by a dent that is one vertical 
position higher than the encoding for the binary bit value 0. Whereas in the 3D encoding each bit 
value is encoded with a fixed two horizontal tile position for encoding, we now have an encoding 
region of size Ik for each bit value. 

As in the 3D case, we utilize a set of decoder tiles (the orange tiles in Figure [4]) to read 
the geometry of the north face of the macro tile to the south. Due to planarity, only one possible 
growth branch is blocked. However, in the case of a bit, the k repetitions of the non-deterministic 
placement of the branching tile makes it highly likely that at least one downward growth will beat 
its competing westward growth over k independent trials, if k is large. In the case of 1 bits, the 
decoder is guaranteed to get the correct answer. 

Therefore, we can analyze the probability that a given bit is correctly decoded by bounding 
the probability of flipping a (biased) coin k times and never getting a single tail. The coin is 
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biased because the south branch of growth must place 3 consecutive tiles before the west branch 
places a single tile to successfully read the bit, which happens with probability 1/^ Thus, by 
making k sufficiently large, we can bound the probability that a given block will make an error by 
incorrectly interpreting a 0-bit as a 1-bit. For a zig-zag system that makes r tile attachments, we 
can therefore set k to be large enough such that with high probability a temperature 1 simulation 
will make all r attachments without error, yielding the following theorem. 

Theorem 2. For any e > and zig-zag tile system T = (T, s, 2) whose terminal assembly has size 
r, there exists a temperature 1, 2D probabilistic tile system that simulates T without error with 

is 



probability at least 1 — e. The scalex of this system is ©(logdj-logr -|- log cr^ log j), the scale 
4, and the tile complexity is 0(|T| logcrr logr -|- |T| log ctt log -). 



Proof. We first note that the tile complexity and scale factor of the simulation are determined 
by the values \T\ and |o"t|; as in the 3D transformation, as well as the choice of the parameter 
k. In particular, the tile complexity of the system is 0{\T\klogaT), the horizontal scale factor is 
scalcx = 0{k\ogaT), and the vertical scale factor is scalcy = 4. Therefore, to show the result we 
need to determine what parameter choice for k is needed to guarantee that all r block placements 
are error free with probability at least 1 — e. 
For a given e > 0, set k equal to 



k = log8/7r-Flog; 



1 



8/7 



As any given block will make an error with probability at most (g)'^, the expected number of total 
errors over r block placements is 



E[errors\ = r{-)^ 



Therefore, for k = logg/^r -|- logg/y ^ the probability of making 1 or more errors is bounded by e, 
and this value of k achieves the theorem's stated bounds. □ 



The north glue coi of the orignal zig-zag tile is 
encoded in the geometry of the output growth on 
the north face of the macro tile. ^ 
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to be 1 . 
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Figure 4: This tile set depicts the macro tiles used to transform a zig-zag system into a probabilistic 
2D assembly system. 



^To make the choice non-biased, multiple copies of the south branching tile types could be included in the tile set, 
making the 3 consecutive placements happen with equal, or even greater probability than the single tile placement 
of the west branch. 
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4 Efficient Assembly of n x n Squares at Temperature 1 



In this section we show how to apply the zig-zag simulation approach from Section [3] to achieve 
tile complexity efficient systems that uniquely assemble nx n squares at temperature r = 1 either 
deterministically in 3D or probabilistically in 2D. 

4.1 Deterministic Temperature 1 assembly of n x n Squares in 3D 

We show that the assembly of an n x n square with a 3D tile system requires at most O(logn) tile 
complexity at temperature 1. 

Lemma 3. For any n, there exists a 2D, temperature r = 2 zig-zag tile system, with north/south 
glueset a of constant hounded size, that uniquely assembles a logn x n rectangle. Further, this 
rectangle can he designed so that a unique, unused glue appears on the east side of the northeast 
placed tile (this allows the completed rectangle to seed a new assembly, as utilized in Theorem^. 

Proof. The binary counter depicted by Rothemund and Winfree^^ is a zig-zag tile set and can 
easily be modified to have a special glue at the needed position. □ 

Theorem 4. For any n>l, there exists a 3D temperature r = 1 tile system with tile complexity 
O(logn) that uniquely assembles an n x n square. 

Proof. The basis of our proof is a construction for a binary counter that is a zig-zag system with 
a 0(1) size north/south glueset. Such a tile system can assemble a logn by n rectangle using 
O(logn) tile types. By Theorem [l| this system can be simulated at temperature 1 in 3D to make 
1 side of an n X n square. By seeding a similar version of this tile system upon the completion 
of the initial assembly, a second and third side of the square can be completed, permitting a final 
filler tile to fill in the body of the square. As the north/south glue set for the counter system is of 
constant size, the same O(logn) tile complexity is achieved for the construction of the square at 
temperature 1 in 3D by combining Lemma [3] with Theorem [T] □ 

4.2 Probabilistic Temperature 1 assembly oi n x n Squares in 2D 

In this section we show that for any value e and positive integer n, there exists a 2D tile system 
with 0(log^ n + logn log ^) tile types that will assemble an n x n square with probability at least 
1 — e. For a fixed constant e, this simplifies to a 0(log^ n) tile complexity. 

Theorem 5. For any e > and positive integer n, there exists a 2D tile system with 0(log^n + 
logn log ^) tile types that will assemble an n x n square with probability at least 1 — e. 

Proof. The basic approach of this theorem is the same as Theorem [4| Consider a zig-zag binary 
counter construction with 0(1) size north/south glueset designed to assemble a logn by n rect- 
angle using O(logn) tile types. We can apply Theorem [2] to this system to get a temperature 
1 probabilistic simulation of the counter. Further, as the counter places at most O(nlogn) tile 
types before terminating. Theorem [2] yields a tile complexity bound of 0(log^n + log n log ^) to 
correctly simulate the counter with probability 1 — e. As with the 3D construction, 3 rectangles 
can be assembled in sequence to build 3 of the borders of the goal n x n square, and finally filled 
with a fill tile type. □ 
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4.3 Kolmogorov Lower Bound for Probabilistic Assembly oi n x n Squares 

Rothemund and Winfree^S showed that the minimum tile complexity for the deterministic self- 
assembly of an n X n square is lower bounded by ^ doj^'iogn ) almost all positive integers n. This 
applies to both 2D and 3D systems. In this section we show that this lower bound also holds for 
probabilistic systems that assemble an n x n square with probability greater than 1/2. 

Theorem 6. For almost all positive integers n, the minimum tile complexity required to assemble 
an n X n square with probability greater than 1/2 is ^{y^y^^)- 

Proof. The Kolmogorov complexity of an integer n with respect to a universal Turing machine U 
is Kjj{n) = min|p| s.t U{p) = hn where hn is the binary representation of n. A straightforward 
application of the pigeonhole principle yields that Kjj{n) > [logn] — A for at least 1 — (^)^ of 
all n (seel^ for results on Kolmogorov complexity). Thus, for any e > 0, Ku{n) > (1 — e) logn = 
O(logn) for almost all n. 

If there exists a fixed size Turing machine that takes as input a tile system and outputs the 
maximum extent (i.e. width or length) of any shape that is terminally produced by the system 
with probability more than 1 /2, then the ^2(log N) bound applies to the number of bits required to 
represent any such tile system that assembles an n x n square with probability greater than 1/2, 
as such a system combined with the fixed size Turing machine would constitute a machine that 
outputs n (note that the unique output of n depends on the probability of assembly exceeding 
1/2 as only one such highly probable assembly can exist in this case). Further, we can represent 
any tileset T = {T, s, r) with 0(|T| log \n\) bits assuming a constant bounded temperature (there 
are at most 4|T| glues in the system, thus all glues can be labeled uniquely with 0(|r|log|r|) 
bits). Thus, for constants d and c we have that dlogn < c|T|log|T|. If we assume T is the 
smallest possible tile set to assemble an n x n square with probability at least 1/2, we know that 
c\T\ log |r| < \T\ log log n as there is a known (deterministic) construction to build any nxn square 
in O(logn) tile types ESI. Thus, we have that dlogn < |r| log logn which yields that the smallest 
number of tile types to assemble a square is |T| = ^{y^^^)- 

The above argument requires the existence of a fixed size Turing machine/algorithm to output 
the dimension of any terminally produced shape assembled with probability greater than 1/2. The 
following algorithm satisfies this criteria. 

OutputDimension( Tile system F = (T, s, r) ) 
• For k from 1 to infinity do the following 

— iterate over all producible fc-tile assemblies of F 

* If the current assembly X is terminal, compute the probability that F assembles 
X. If the probability exceeds 1/2, output the dimension of X and halt. 

□ 

4.4 Why not 0{ ^^°f^^ ) square building systems at Temperature r = 1? 

Temperature r = 2 assembly systems are able to assemble nxn squares in 0{ ^^°f^^ ) tile complex- 
ity, which meets a corresponding lower bound from Kolmogorov complexity that holds for almost 
all values n^. Further, this can be achieved with zig-zag tile systems. So, why can't we get 
^( lo'giogra ) complexity for temperature r = 1 assembly in 3D? One answer is, we don't know 
that we can't. However, the direct application of the zig-zag transformation to the temperature 
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2, 0( 

log log n ) ^^^^ complexity system does not work. In particular, the temperature r — 2 result is 
achieved by picking an optimal base for a counter, rather than counting in binary. In general, a 
k X n rectangle can be assembled by a zig-zag tile system in 0{k + ) tile complexity by imple- 
menting a k digit, base nfc counter IS. However, inherent in counting base n^, the counter uses at 
least nk north/south glue types. Thus, the zig-zag transformation yields a tile complexity of 

(/c + n fe ) log(n ) = r2(log n) 

Therefore, regardless of what base we choose for the counter, we get at least a Q{logn) tile 
complexity. The same problem exists for the base conversion scheme used in E . It is an interesting 
open question whether or not it is possible to achieve the optimal O( io'gk)gn ) complexity bound 
at temperature r = 1. 

5 Turing Machine Simulation for General Computation 

By extending non zig-zag Turing machines from the literature we show how to implement a 
zig-zag tile set to simulate the behavior of an arbitrary input Turing machine. From Theorem [T} 
we can simulate this system at temperature 1. The details of the construction are given in the 
Appendix in Section [Bj 

In this section we discuss our results for temperature 1 tile assembly systems capable of sim- 
ulating Turing machines. Our results are based on existing temperature 2 self-assembly systems 
capable of simulating Turing machines ESI. With straightforward modifications, these construc- 
tions can be modified into zig-zag systems, and thus can be simulated at temperature 1 in 3D 
or in 2D probabilistically. Our first Lemma states that such temperature 2 zig-zag systems exist. 
For simplicity of bounds, we assume the Turing machine to be simulated is of a constant bounded 
size (constant bound number of states, alphabet size) and refer the reader to the appendix for tile 
complexities in terms of these variables. 

Lemma 7. For a given Turing machine T there exists a temperature 2 zig-zag tile assembly 
system that simulates T. More precisely, there exists a zig-zag system that, given a seed assembly 
consisting of a horizontal line with north glues denoting an input string, the assembly will place a 
unique accept tile type if and only if T accepts the input string, and a unique reject tile type if and 
only ifT rejects the input string. 

By combining this Lemma with the zig-zag simulation lemma, we get the following results. 

Theorem 8. For a given Turing machine T there exists a temperature 1, 3D tile assembly system 
that simulates T, even when assembly is limited to one step into the third dimension. Thus, 3D 
temperature 1 self-assembly is universal. 

Proof. This theorem follows from Lemma [7] and Theorem [l} The exact tile complexity and scale 
of the simulation in terms of the number of states and alphabet size of the Turing machine are 
derived in the Appendix in Section [B} □ 

Theorem 9. For a given Turing machine T there exists a temperature 1, 2D probabilistic tile 
assembly system that simulates T. In particular, for any e > and machine T that halts after N 
steps, there exists a temperature 1, 2D system with 0(log + log ^) tile types in terms of N and 
e that simulates T and guarantees correctness with probability at least 1 — e. 
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This theorem follows from combining Lemma [7] and the construction for zig-zag simulation 
used in Theorem [2] with parameter k = log As discussed in Section 3.4 the probability of an 



single error in simulation of a given zig-zag set is 1/2 . Thus, the expected number of errors for 
k = log Y becomes implying that the probability of 1 ore more errors is bounded by at most e. 

The exact tile complexity and scale of the simulation in terms of the number of states and 
alphabet size of the Turing machine are derived in the Appendix in Section [B} 



6 Simulation Results 

We have developed automatic conversion software to generate temperature 1 3D or 2D systems 
given any input zig-zag tile set. We have further verified the correctness of our algorithms and 
converter software by implementing a number of tile systems such as Sierpinski tile sets and binary 
counters. We have run the simulations on Caltech's Xgrow simulator and ours, the system 
performs well. 




Figure 5: Snapshots of Sierpinski in 3D Zig-Zag Tile System 

We are interested in experimentally comparing the fault tolerance of our temperature 1 prob- 
abilistic system with other fault tolerance schemes such as snaked proofreading systems. We have 
implemented both systems with similar tile complexity constraints to see under what settings ei- 
ther system would begin to make errors. We have done preliminary simulations of the 2D systems 
in the kTAM with Caltech's Xgrow. Our test case was a simple parity checking tile system and 
our metric was whether or not the correct parity was computed. Our preliminary tests used the 
input string " 1000" . The block size of the snaked system is 6 x 6, and the length of each hit{K) of 
probabilistic assembly tiles is 5. The glue strength of the probabilistic assembly tiles are changed 
to 2. 

We used Gmc = 13.6, Gse = 7.0. Both tile systems outputted correct results. But as we 
increased the Gse, the performances of the two systems began to differ. The snaked system became 
unstable, facet errors occurred and the tile mismatches increased. In contrast, the probabilistic 
assembly tile set grew faster and the output of the system was still correct. We used the r = jf^ 
ranging from 2 to 0.5, and the results show that the 2D probabilistic assembly system is more 
stable than the snaked system within the parameters we tested. 
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The comparison of the probabihstic and snaked systems is very prehminary because the two 
system use different temperatures. We adjusted the glue strength of the probabihstic system to 
make it comparable with the snaked system. To more accurately compare the two systems, we 
need to investigate other methods and test sets, both for probabilistic and snaked systems. 

7 Further Research Directions 

There are a number of further directions relating to the work in this paper. The most glaringly 
open problems are related to the power of deterministic temperature 1 self-assembly in 2D. For 
example, it is conjectured that the optimal tile complexity for the assembly of an n x n square 
at temperature 1 in 2D is 2n — 1 tile types. While it is straightforward to achieve this value, the 
highest lower bound known is no better than the small 0{ ^^°y^^ ) lower bound from Kolmogorov 
complexity. Similarly, little is known about the computational power of this model. It seems that 
it cannot perform sophisticated computation, but no one has yet been able to prove this. 

Another direction is the concern that temperature 1 systems are prone to multiple nucleation 
errors in which many different growths of tiles spontaneously attach without requiring a seed. One 
potential solution is to design growths such that assemblies that grow apart from the seed have 
a high probability of self-destructing by yielding growths that block all further growth paths for 
the assembly. In contrast, properly seeded growths by design could have the self-destruct branches 
blocked by the seed base they are attached to. 
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A A Simple 3D Example 

See Figure [6] and [7} 
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Figure 6: The tileset in (a) is a zig-zag tileset that swaps the right most a to b every other hne. 
In (b), (c), and (d), a cohection of macro tiles are provided, according to Theorem [TJ This set will 
simulate the given tileset system using 3 dimensions at temperature r = 1. 
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B Constructing a Turing Machine with Zig-Zag Tile Sets at r = 2 
in 2D 



B.l Converting 

A Turing machine is a 7-tupleEa, {Q,^,T,6 

: QOi Qaccept: Qreject) : where Q, S, F are all finite sets 

and 

1. Q is the set of states, 

2. E is the input alphabet not containing the blank symbol ^, 

3. r is the tape alphabet, where ^ G F and S C F, 

4. 5 : Q X T — > Q x T x {L, R} is the transition function, 

5. qo £ Q is the start state, 

6. Qaccept G Q is the accept state, and 

7. Qreject £ Q is the reject state, where Qaccept 7^ Qreject- 

All of the glues at the side north or south of the tile regard as the alphabet in set T. The gn of 
tile is the output and gs is the input in our turing machine, so the action of placing one tile with 
difference gn and gs on the supertile regard as one write operation, placing one tile with the same 
gn and gs regard as one read operation. 

Basically, the tiles on the west of growth direction will be used to the main calculation, and 
the tiles on the direcition of east will be used to copy the glues from south to north. 

First, let's consider the equation S{q,c) = {q',c',R). The head of TM is located at the c 
position and the current state is q. Then the record of the current position of tape is changed to 
c', the current state is changed to q' , and head move to right (R). Three types of the tiles, which 
are indicated by to, ti, and t2 in Figure [8|3, are used to simulate this process. The glue notation 
at the south of the tile to means that the head is at the column of tile to and input alphabet is c. 
The glue notation (q) at the west of the tile to means that the current state is q. The glues at the 
north and east indicate the output and next state seperately. The east glue {q'{R}) of the tile to 
indicates the next state {q') and earring the extra information to notify the next tile that the head 
is coming. 

The equation 5{q, c) = {q' , c' , L) indicate that the head will move to left. Figure [8^ shows how 
to move the head to left. Again, the head is at the column of tile to, the current state is q, and 
input is c. Then the tile to output north glue c'{L}, and east glue g'. The tile ti will translate 
the glue information (c'{L}) to state information [q'r{L}) which contains the command of moving 
left. At last, the tile t2 add the head notation to the glue of the current column. 

If the head don't move after read one character from tape, the corresponding tile is to showed 
in Figure [S}:. Figure [8]i and Figure show the tiles at the left and right sides of zig-zag structure 
seperatly. Figure |8f shows the seed bar, which initializes the turing machine, includes the start 
state go; the input sequence (a, 6, ... , z) in the type etc. 

B.2 Complexity Analysis 

In the case of r = 2 2D zig-zag tile system, the number of the tile types for state transfer (tg in the 
Figure p| is \5\. And the number of the auxiliary tile types for state transfer (ti and t2 in Figure 
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Figure 8: Zig-Zag Tile Structure for Contructing Turing Machine. (a) 5{q, c) = 
(g',c',L);(b) 6{q,c) = {q',c',R); (c) 6{q,c) = (g',c',0); (d) Left side; (e) Right side; (f) Seed bar: 
go and input sequence a,b, . . . ,z 



[8)3 ) is \T,\ X IQI X 2 (move left) and x |Q| x 2 (move right). The number of tile ti in Figure 
does not need be included because it is counted in the auxiliary tile types {t2 in Figure [sja ) . 
The other tile types are the tiles for transfer alphabet information, the number of it is |(5| x |S| x 2 
(2 indicate the east direction and west direction of the growth). 

There's also need some special types of tile in both sides of the zig-zag structure to change the 
growth direction, it will cost 4|Q| tiles. If the length of the tape is n, then the number of the tile 
types for constructing seed bar will cost n + 2 tiles (See Figure [8]F). 

So the total tile types is T = \6\ + |S| x \Q\ x 2 + |S| x \Q\ x 2 + |Q| x |S| x 2 + 4|Q| + n + 2 = 
|<5| +6|Q||S| + 4|Q| + n + 2. The upper bound of tile types is 0{\6\ + +n). Considering that 

IQI < \S\ < IQPiSp and n < the complexity of tile types is 0(|<5| + |Q||S| +n) = 0{\6\). 

Each of the state transfer cost two lines of the zig-zag structure (east direction line and west 
direction line above). Given a input tape and start state, if the number of the state transfers to 
be passed before it stop is r, then the lines of the tile structure will be 2r. So the space used in 
calculation is 0{nr). 

In the case of r = 1 3D zig-zag tile system, the length of each mapped tile depend on the length 
of encoding code of glues. If the number of glues of the 2D zig-zag tiles is G, then G = 0{\6\). 
The complexity of tile types in 3D is 0{\6 \ log G) = 0(|5| log |5|), the complexity of space is 
0(nr log G) = 0{nr log \5\). 
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In the case of r = 1 2D probabilistic zig-zag tile system, the length of each mapped tile depend 
on the parameter K and the length encoding code of glues. The parameter K of the r = 1 2D 
probabilistic zig-zag tile system, is a constant depend on the probability of correction of result. 
K = 0(logr + logi). The space complexity is 0{nrK\o^G) = 0{nr\og\5\), and the complexity 
of tile types is 0(|(5| log |(5|). 



Table 2: The Complexity of Zig-Zag Turing Machine 





r = 2 2D 


T = 1 3D 


r = 1 2D Prob. 


Space 


0{nr) 


0{nr log \5\) 


0{nr log \6\) 


Tiles 


oi\s\) 


0{\6\log\6\) 


0{\6\log\6\) 



\6\: the number of state functions; 
n: the length of the tape (number of alphabet); 
r: the number of the state transfers to be passed before it stop; 
K: the parameter of probabilistic zig-zag tile system; 

B. 3 Algorithms 

The detailed algorithms for transfering any turing machine to zig-zag tile set at temperature r = 2 
in 2D is listed in Algorithm [l| 

C Convert Arbitrary Zig-Zag Tile Set from r = 2 to r = 1 

C. l Notation 

A tile t is a four sided Wang tile denoted by the quadruple (gn, Qe, gs, Qw), g-n, Qe, gs, 9w denote the 
glue type of the four sides: North, East, South, and West. 

C.2 The Converting Table 

When the Zig-Zag tiles in 2D are categorized into some types, the tile types can be mapped to 3D 
and 2D probabilistic tile sets directly. Table [3] list all of the relationships between the 2D tiles and 
3D tile sets. 

The tiles (see the 3D figures in Table [s]) are different from each other. The lines between two 
adjacent tiles denote the glues which strength are 1 and they are also different from each other 
except adjacent glues. Large squares denote the tiles in the plane z = 0, and small squares denote 
the tiles in the plane z = 1. The encoded code(en) showed in the figures in Table [3] is two bits 
long(maxbits=2). The parameter K of Probabilistic Zig-Zag is 2 in Table [3j 
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Algorithm 1: Turing-machine-to-zigzag(): Convert TM to Zig-Zag Tile Set at r = 2 in 2D 
Input: go. the start state 
str, the input string at the tape 
5, the state transition functions 

Output: T2d,2t, zig-zag tile set at temperature 2 in 2D 

Q (() ; /* All of the states of the 5. */ 

Yj ^ <j) \ /* All of the input/output of the 5 and the alphabet at the tape. */ 

foreach / G 5 do 

Q <— Q U { current state of /} ; 
Q ■(- Q U { next state of /} ; 
S<(-EU{ input of/} ; 
S <(- S U { output of /} ; 

foreach c € str do 

L E ^ S U {c} ; 

foreach g e Q do 

/* A four sided Wang tile denoted by the ordered quadruple 

(north(t ) , east {i), south(t ) , west {t)) */ 
T2d,2t ^ T2d,2t U {fAle{.sl{q},qr, sltr, (j))} 
T2d,2t ^ T2d,2t U {tile{sltr, q, sl{q}, (f))} ; 



f2d,2t 



U {tile{sr{q}, (j), srtl, q)} ; 



/* left, down */ 

/* left, up */ 
/ * right , down */ 
/* right , up */ 



T2d,2t 

T2d,2t ^ T2d,2t U {tile{srtl, (f>, sr{q}, qr)} 
foreach c G S do 

T2d,2t ^ T2d,2tU {tile{c,q,c,q)} ; /* tiles for coping */ 

T2d,2t ^ T2d,2t^ {tile{c,qr,c,qr)} ; /* tiles for coping, return to left sides. */ 
T2d,2t ^ T2d,2t^ {tile{c{H},qr,c{H},qr)} ; /* Right moving auxiliary tiles: ^2 */ 
T2d,2t ^ T2d,2t^ {tile{c{H},q,c,q{R})} ; /* Right moving auxiliary tiles: ti */ 
T2d,2t T2d,2t^ {tile.{c{H},qr{L},c,qr)} \ /* Left moving auxiliary tiles: t2 */ 

foreach / € 5 do 

switch head moving do 
case LEFT 

T2d,2t ^ T2d,2t U {tile{f.output{L}, f.stateout, f-input{H}, f.statein)} ; 
T2d,2t T2d,2t U {tile{f .output, f.stateoutr, f.input{L}, f.stateoutr{L})} ; 

case RIGHT 

|_ T2d,2t ^ T2d,2t U {tile{f. output, f.stateout{R}, f-input{H}, f.statein)} I 
case No moving 

|_ T2d,2t T2d,2t U {tile{f.output{H}, f.stateout, f.input{H}, f.statein)} ; 



/* Seed Bar 

T2d,2t ^ T2d,2t U {tile{srtl, cf), cp, gr)} ; 
T2d.2t ^ T2d,2t U {tile{sl{qo},gl,<j),(l))} ; 
for i <— to {\str\ — 1) do 
c <— ith alphabet of str ; 
if z = then 
L T2d,2t ^ T2d,2t U {tile{c{H},Si,^,gl)} ; 

else if i = {\str\ — 1) then 
L T2d,2t ^ T2d,2t U {tile{c, gr, (j), Sj-i)} ; 

else 

L T2d,2t ^ T2d,2t U {tile{c, Si, (j), Sj_i)} ; 
return T2d,2t; 



/* to */ 
/* ti */ 

/* to */ 
/* to */ 

*/ 

/* SEED */ 
/* SEED_L */ 
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Table 3: Zig-Zag Tile Set Mapping 



Zig-Zag in 2D T = 2 Zig-Zag in 3D r = 1 Prob. Zig-Zag in 2D r = 1 



Direction West, West 1 
(DWWl) 



Encoding Area 



,\ 1 Encoding Area M5B 



- (y,-) 



|^gjn = (x,a) 



Direction East, East 1 
(DEEl) 



MSB Encoding ^ 



II ^1 



ajn . II, a) 



Encoding Area 



IT 

E 



Turn at East, West 1 
(TEWl) 




E 



Turn at West, East 1 
(TWEl) 



f^SB Encoding 




MSB Encoding Area 



S3 



Continued on next page 
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Zig-Zag in 2D T = 2 Zig-Zag in 3D r = 1 



Prob. Zig-Zag in 2D r = 1 



B 



Direction East, North 2 
(DEN2) 



B 



Direction West, North 2 
(DWN2) 



|«.gjn.tt. 



B 



Side West, North 2 (SWN2) 



vm 



B 



Side East, North 2 (SEN2) 



S13 



Encoding Area MSB 



:S0 



Fixed Location, Direction 
West (FW) 



"1 



--l-> 






;-> 


» -in — 


- ♦ 1 




-— 1 


l ° ' 


LSB MSB Enco 



LS3 



i 



SO: 



MSB 



:si 



Fixed Location, Direction s„ 
East (FE) 
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Continued on next page . . . 
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Zig-Zag in 2D T = 2 Zig-Zag in 3D r = 1 



Prob. Zig-Zag in 2D r = 1 



^ Direction 
West, West 2 (DWW2) 



T_0 








I Encoding Area MSB LSB 
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1 
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Turn at West, West 2 
(TWW2) 
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Continued on next page . . . 
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Zig-Zag in 2D T = 2 Zig-Zag in 3D r = 1 



Prob. Zig-Zag in 2D r = 1 



MSB Encoding Area 



IT 

E 



Turn at West, East 2 
(TWE2) 



Hi l» 



< 




< 




















— H 









— -1 


1 — 1 


'■-h 


I—- 
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C.3 Algorithms for Converting Tile Type 

To converting arbitrary zig-zag tile set, the first step is to categorize the tiles into the sixteen types 
listed in Table [3j Then the zig-zag tiles in 2D can be converted directly to zig-zag tile set in 3D or 
probabilistic zig-zag tile set in 2D. We can recgonize all of the tile types in 2D by observing the 
positions the tiles be placed during the growth. In an other words, a simulator have to be used to 
get the types of tiles in 2D. 

The tiles are categorized into sixteen types, it seems it will cost much effort to implement the 
converter by software. As we apply some tricks on the algorithms, the creating of tile sets converter 
would be simple. The tile set type DWWl, DWW2, TEWl, TEW2, TWW2 have similar structure, 
the differences between those types are the positions of the input and output glues. While the tile 
set type DEEl, DEE2, TWEl, TWE2, TEE2 also have the similar structure, they are the mirror 
of the tile set type DWWl, DWW2, TEWl, TEW2, TWW2. And SWN2 is the mirror of SEN2, 
DWN2 is mirror of DEN2. 



C.3.1 Categorize Zig-Zag Tiles 

A simulator is used to detect the tile types in this algorithm. The only restriction is that all of 
the tile types should be occured at least once in resonable steps during simulating. The algorithm 
will return right after all of the tile types are detected. The algorithm is listed in Algorithm [2] 



C.3. 2 Zig-Zag in 3D 

The algorithms to create the tile sets in 3D are listed in Algorithm |3j|4l[5}[6K7K8K9Kl0KTl}[T2 



13 , 14 The number of the binary bits used in the codes is denoted by maxbits. 



C.3. 3 Probabilistic Zig-Zag in 2D 

The algorithms for converting the zig-zag from temperature r = 2 to temperature r = 1 are similar 
to that in 3D. Using the same algorithms to encode all of the glues at the north or south of tiles 
with strength 1 (See Algorithm [3| |4]) . 

The decoding tile sets are different from the zig-zag in 3D. The parameter K is introduced in 



the probabilistic zig-zag tile set. Figure 11, 12 shows the tile set for detecting one bit of the code 
by using K{= 4) groups of the detect tile set. Figure [l3] shows one of the complete decoding tile 
sets which have similar function as showed in Figure [9j 

The Encoding Area is a bit different from that in 3D. The length for each bit of the code 
in the Encoding Area is depend on the parameter K. The length of the Encoding Area will be 



28 



Algorithm 2: Categorize-Zig-Zag-tiles() 



Input: T2^2f Un-categorized tile set at temperature 2 in 2D 
Output: T2d,2t< Categorized tile set at temperature 2 in 2D 

numTypes <— ; 

while numTypes < |r2d2tl 

select one of the tile t' G 2t '^^^^^ attach to the supertile at position pos ; 
record this step {t' , pos) and append to array steps ; 
if the type oft' is unknown then 
adjidx i- ; 

numTypes numTypes + 1 ; 

for i ^ {\steps\ — 1) to 1 do 

if pos is adjacent to steps[i].pos and the glue between them are the same then 
adjidx i ; 
adjTile ^ steps[i].t ; 

dir the side of t' that is adjacent to steps[i].t; 
Goto end of this for loop ; 
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if adjidx is between {\steps\ — 1) and 1 then 
switch the value of dir do 
case EAST 

if strength oft'-g^ > 1 then 

if strength oft'.gn > 1 then The type oft' <(- SWN2; T2d,2t ^ T2d,2t U{t'} ; 
else if strength oft'.gn, > 1 then The type of t' <- FW; T2d 2t <- 2t 
U{f} ; 
else 

if strength oft'.gn > 1 then The type of t' ^ DWN2; T2d.2t ^ 72<i,2t U{t'} ; 
else if strength of t' .g^ > 1 then The type of t' ^ DWW2; T2d 2t ^ T^d 2t 
U{t'} ; 

[_ else The type of i' <(- DWWl; T2d,2t ^ 72d,2t U{i'} ; 

case SOUTH 

switch the type of adjTile do 
case DEN2 or SEN2 

if strength oft'.gg > 1 then The type of t' TEE2; T2d2t ^ T2d2t 
iJ{t'} ; 

else if strength oftf.g^ > 1 then The type oft' TEW2; T2d,2t 
T2d.2t U{t'} ; 

else The type of t' ^ TEWl; T2d,2t ^ T2d,2t U{t'} ; 
case DWN2 or SWN2 

if strength of t'.ge > 1 then The type of t' TWE2; T2d 2t ^ T2d 2t 

U{t'} ; 

else if strength oft'.g^ > 1 then The type of t' <(- TWW2; T2d,2t ^ 
T2d.2t U{i'} ; 

else The type of t' ^ TWEl; T2d,2t ^ T2d,2t U{t'} ; 

case WEST 

if strength oft'.g^ > 1 then 
if strength oft'.gn > 1 then 

I The type of t' ^ SEN2; T2d,2t T2d,2t U{t'} ; 
else 

L The type of t' ^ FE; T2d,2t ^ T2d,2t U{t'} ; 

else 

if strength oft'.gn > 1 then The type of t' ^ DEN2; T2d,2t ^ T2d,2t U{t'} ; 
else if strength of t'.ge > 1 then The type of t' DEE2; T2d 2t ^ T2d 2t 
U{t'} ; 

[_ else The type of t' ^ DSBL; T2d,2t ^ 72d,2t U{t'} ; 
otherwise Error, Ignored ; 



return T2d,2t: 



Algorithm 3: Tileset-convert- all-categories () 
Input: T2d.2t, Categorized tile set at temperature 2 in 2D 
Output: T^d,it, Tile set at temperature 1 in 3D 
/* 

T2d,2t — Tdwwi U TdwW2 U Ttewi U TteW2 U TtwW2 
U TuEEl U TeiEE2 U Ttwei U TtwE2 U TteE2 
U TpB U TfIV U T5VKiV2 U TsEN2 U TdwN2 U TdeN2 

*/ 

Tsd.it ^ 4> \ /* The tile set in 3D temperature 1 */ 

Gns ■^01 /* The glues at the north and south sides of the tile */ 

S2dat ^ {si\si = Strength of all of the glues of tj, tj e T2d,2t } ; 
foreach U e T2d,2t do 
\i 1 = Sti.g^ then 

L ^ns ^ Gns U {ti.gn}\ 

if 1 = Sf-.g^ then 
|_ Gns ^ Gns U {U.gs}; 

Encode the glues in set Gns by binary codes ei\i E Gns', 

£ ^ G Gns} /* £ contains all of the code of glue £ Gns- */ 

maxbits ^ [log|G„5|]; 

Tsd.lt ^ T3d,lt U GENERATE-DECODE-TILE-SET(T<i„ie/t UTd„„3ht,S'2d,2t,f ■ maxblts) ; 
TsdM ^ Tsrf.lt U GENERATE-ENCODE-TILE-SET(rdjrie/t U Tdirright, 5*2^,24, maxbits) ; 

return T:>,d,u\ 



2K X maxbits. The mapping between the zig-zag and probabihstic zig-zag for each of the tile types 
is showed in Table [3j the implemention use the similar algorithms as that used in 3D. 

The success ratio of constructing zig-zag structure depends on the parameter K, but we noticed 
that it also depends on the number of zero in encoding code, because there exist false positive in 
detecting the zero bits of the encoding codes. We can select the codes which contain many one 
bits for the encoding code to reduce the error ratio. 

The algorithms for probabilistic zig-zag in 2D are omited. 
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Algorithm 4: Generate-decode-tile-set() 



Input: T2d,2f, Categorized tile set at temperature 2 in 2D 
S2d,2t, all of the strength of glue In T2d,2t 
£, all of the code of glue G Gns 

maxbits, the number of binary bits to encode all of the glues G Gns 
Output: Tg^ Tile set at temperature 1 in 3D 

G2w ^ </> ; /* The glue set of input from east to west */ 

G2e ^ <f> \ /* The glue set of input from west to east */ 

foreach U G Tdwwi U Tdww2 U Tuwn2 do 

L G2w ^ G2w U {U.ge}; 
foreach ti e Tdeei U Tdee2 U Tden2 do 

foreach € G2w do 

foreach e Tdwwi U T'£)tyvy2 U Tdv^jv2 do 
if tj .ge = gi and St^ .g^ < 2 then 

/* St-.g, is the strength of g^ of tile tj, St-.g^ G S2d,2t */ 

Gw,i ^ Gw,i U {tj.Qs} 

It = It ^ Generate-decode- TILE-TO- WEST(5j, G^,,, f , maxbits); 
foreach e G2e do 

Ge,i = (/>; 

foreach tj e Tdeei U Tde£:2 U Tden2 do 
if fj..9ui = gi and stj.g^ < 2 then 

/* St,.g, is the strength of g^ of tile tj , St^,g^ € S'2d,2t */ 

Ge,i Ge,i U {tj.gs} 
^3d,lt ^ ^3d,lt 'J GENERATE-DECODE-TILE-TO-EAST(5i,Ge,i,^,maa;6iis); 

return T^^ -^^; 
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Algorithm 5: Generate-decode-tile-to-west() 



Input: gin, the input glue of the current tile set 

Gw,i, the glue set to be the output glues of current tile set 

£, all of the code of glue e Gns 

maxbits, the number of binary bits to encode all of the glues e Gns 
Output; Tg^ j^j, Tile set at temperature 1 in 3D 

/* Contruct a binary tree according to the encoding code of each it ems in G^.i • */ 
for each e G^a do 
curnode root; 

foreach hit of Cg^ from MSB to LSB do 
if bit = 1 then 

if curnode have no right child then 
[_ create right child of the curnode; 

|_ curnode ^ curnode. right_child; 

else 

if curnode have no left child then 
[_ create left child of the curnode; 
L curnode ^ curnode. Ieft_child; 

Traversal the tree by pre-order algorithm: Part of the tiles set will be generated and saved to Tg^ in 
each visitation. The input glue of the tile set is g^. See Figure[9] 
return TL, , , ; 



Algorithm 6: Generate-decode-tile-to-east() 



*/ 



Input: gin, the input glue of the current tile set 
Ge,i, the glue set to be the output glues of current tile set 
£, all of the code of glue e Gns 

maxbits, the number of binary bits to encode all of the glues e Gns 
Output; Tg^ j^j, Tile set at temperature 1 in 3D 

/* Contruct a binary tree according to the encoding of each items in Ge^i 
foreach e G^j do 
curnode ^ root; 

foreach bit of Cg. from MSB to LSB do 
if bit = 1 then 

if curnode have no right child then 
[_ create right child of the curnode; 

curnode ^ curnode. right_child; 
else 

if curnode have no left child then 

[_ create left child of the curnode; 
curnode -s— curnode. Ieft_child; 

Traversal the tree by pre-order algorithm: Part of the tiles set will be generated and saved to 



each visitation. The input glue of the tile set is gin- See Figure 10 



3d,lt 



in 



return T'; 
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Algorithm 7: Generate-encode-tile-set() 



Input: T2d.2t, Categorized tile set at temperature 2 in 2D 
£, all of the code of glue G Gns 

maxbits, the number of binary bits to encode all of the glues e Gns 
Output: Tg^ Tile set at temperature 1 in 3D 

UdAt = 

It ^ It ^ GENERATE-ENCODE-TILE-SET-TO-EAST(T2d,2t, maxUts); 

Tg^ It ^ T'sd It U GENERATE-ENCODE-TILE-SET-TO-WEST(T2d,2t, maxbits); 

^3d,lt ^ ^3d,lt GENERATE-ENCODE-TILE-SET-OTHERS(T2(i,2t,^,'Tiaa;&iis); 

return TL. 



Algorithm 8: Generate-encode-tile-set-to-west() 
Input: T'2d,2t. Categorized tile set at temperature 2 in 2D 
£, all of the code of glue e Gna 

maxbits, the number of binary bits to encode all of the glues € Gns 
Output: Tg^ Tile set at temperature 1 in 3D 

^3d,it = ^; 

foreach U &Tdwwi 

9in ^ {U-ge,ti.gs)\ 
Qout ^ {ti-9w, —)', 

T^d,it ^ ^M.it U GENERATE-ENCODE-TiLE-TO-WEST(DW^Wl,5i„, gout, etj.g^,maa;6its) ; 
/* C't^.gr, is the encoding code of glue Qn of tile ti */ 
foreach ti S Ttewi do 

9in ^ ^i'Qst 

gout {ti-gw, 

Tg^ If ^ T^d It U GENERATE-ENCODE-TiLE-TO-WEST(riJ14^1,5ij„,5o„t,ef;.g„,maa;tos); 
foreach ti G Tdww2 do 

5m {ti-ge,ti-gs)\ 
gout ^ -gw T 

^3d,u ^ ^3d,it '-^ Generate- ENCODE- TILE-TO- WEST(£>l^W2,gi„,5o„t,eti.g„,max6its); 
foreach ti G Ttww2 do 

gin ^2 -^s I 

gout ii-gwi 

T^d,lt ^ ^3d,lt GENERATE-ENCODE-TILE-TO-WEST(TWW2,5i„,5fo„t,et..g^,maa;6its); 

foreach e Ttbh'2 do 

5in ^ U-gs: 
gout ^i^gwi 

Tg^ It ^3d It ^ Generate-encode- TILE-TO- west(T£;VF2, gin, gout, eti.g„ , maxbits); 
return Ti. , + ; 
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Algorithm 9: Generate-encode-tile-set-to-east() 

Input: T2d,2t, Categorized tile set at temperature 2 in 2D 
£, all of the code of glue G Gns 

maxbits, the number of binary bits to encode all of the glues € Gns 
Output: Tg^ Tile set at temperature 1 in 3D 

foreach ti £ Toeei do 

Qin ^ {^i -dw 7 -Ss ) I 
Qout ^ {ti-ge, — )\ 

^3dU'^^3dlt^ GENERATE-ENCODE-TILE-TO-EAST(Z)i;£l,fifj„,5o„t,e(^.g„,maa;6iis) ; /* eti.g, 

is the encoding code of glue Qn of tile */ 
foreach ti e Ttwei do 

Sout {ti-9e, 

^3d,it ^ ^3d,it ^ Generate-encodE'TILE-to-east{TW El, Qin, gout, eti.gr,, maxbits); 
foreach U &Tdee2 do 

gin ^ (ti-gwi^i-gs) ' 
gout ti.ge', 

Tg^ It ^ It U Generate-encode-tile-to-east(£'£;£;2, gi„,5io«t, et^.g^,, maxbits); 
foreach e Ttee2 do 

,9ot/,t i'i-.9el 

^3d,lt ^ ^Sd.lt 'J GENERATE-ENCODE-TILE-TO-EAST(T£;i;2,5j„, flout, et4.Sn>"^'^2;6jts); 

foreach ti G Ttwe2 do 

gin -^/s » 

gout ^ ii-5el 

Tg^ It ^ "Tad It ^ Generate-encode-tile-to-east(TW£'2, gin, gout, ^ti.g,, , maxbits) ; 
return Ti. .^; 
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Algorithm 10: Generate-encode-tile-set-others() 



Input: T2d,2t, Categorized tile set at temperature 2 in 2D 
£, all of the code of glue G Gns 

maxbits, the number of binary bits to encode all of the glues e G„s 
Output: Tg^ Tile set at temperature 1 in 3D 

foreach ti e Tuwn2 do 
Qin ^ {U-ge,ti-gs)\ 

Qout ^ ^i-9nf 

"^sdAt ^ "^sdAt U Generate-connectioi<i-tile{DW N2, Qin, Qout, maxbits); 
foreach ti G Tden2 do 

9in ^ ^i'Qwi^i'Qs) y 

Tg^ It ^ T^d It U GENERATE-CONNECTiON-TiLE(£'ii;Af2, 5j„, go„t, maxbits); 
foreach ti G T'siyjv2 do 

Qin ^ ti.ge, 
gout -gn f 

Tg^ It ^ It U Generate-connection-tile(S'VFA/'2, gin, gout, maxbits); 
foreach ti G Tsen2 do 

i?out ii-g-ni 

Tg^ If ^srf It U GENERATE-coNNECTiON-TiLE(S'i;A''2, ^-j^, gout, maxbits); 
foreach ti G Tpy\r do 

gin ^i-gei 

gout ^ ii^gwf 

T^d,it ^ T^d,it^ Generate-encode-tile-fixed(FW, 5(j„, gout, eti.g„, maa;6its) ; /* et,.g„ is 
the encoding code of glue Qn of tile ti */ 
foreach ti &Tfe do 

gin ^ ^i^gwi 
gout ^ ti-ge'< 

Tg^ It ^ It GENERATE-ENCODE-TILE-FIXED(i^i5,5l„,5o„t,efi.g„,maa;tos); 

return 
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Algorithm 11: Generate-encode-tile-to-west() 



Input: tiletype, the type of the tile 

gin, the input glue 

gout, the output glue 

e„, the code of the north glue 

maxbits, the number of binary bits to encode all of the glues G G„s 
Output: Tg^ j^j, Tile set at temperature 1 in 3D 

/* Generate distinct tiles as showed in Table [sj. */ 
if tiletype = DWWl then 

Generate the tiles t'^, t[,t'2,t'^, t'^, tg and put it to Tg^ j^^; 

/* The ^Q, t']^, ^3, ^4, ^5 are denoted separately by t_0, t_l, t_2, t_3, t_4, and t_5 
in the figure of DWWl in the Table [3]. The glue is go^t at the west of the 
tile tp. */ 

else if tiletype = DWW2 then 
1^ Generate the tiles ip, t'l and put it to Tg^ ^f.; 

else if tiletype = TEWl then 
1^ Generate the tiles t'Q,t[,t'2,t'^,t'^,t'^ and put it to T^^^^; 

else if tiletype = TEW2 or tiletype = TWW2 then 
1^ Generate the tiles if,, t[ and put it to Tg^ 

Generate the tiles in the Encoding Area and put it to Tg^ 

/* The positions of the tiles are depend on the e„ ; Each bit of the e„ are encoded 
by two tiles in the plauie z = l; The tiles will place at the position 'l'(the 
dotted line denoted by '1' in the figures if the bit is 1, while the tiles will 
place at position '0' if the bit is 0; The encoding of the most significsint bit 
(MSB) of the e„ is place at the left side of the Encoding Area, and the least 
siginif leant bit(LSB) of the e„ is place at the right side of the Encoding Area. 
All of the encoded tiles are connected by the tiles in the plane z — 0. */ 

if tiletype = DWWl or tiletype = DWW2 or tiletype = TWW2 then 
1^ Generate the tiles between glue A and glue gin, put those tiles to Tg^ 

else if tiletype = TEWl or tiletype = TEW2 then 
1^ Generate the tile with glue g^n, put it to Tg^ 

return T^^^^,; 
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Algorithm 12: Generate-encode-tile-to-east() 



Input: tiletype, the type of the tile 

gin, the input glue 

gout, the output glue 

e„, the code of the north glue 

maxbits, the number of binary bits to encode all of the glues e Gns 
Output: Tg^ j^j, Tile set at temperature 1 in 3D 

if tiletype = DEEl then 

Generate the tiles t'Q,t'i,t2,t'^,t'^,t'^ and put it to Tg^ 

/* The <Q, t'j^, tj, tg, ^5 are denoted separately by t_0, t_l, t_2, t_3, t_4, and t_5 
in the figures of Table [sj. The glue is gout at the east of the tile t'^. */ 

else if tiletype — DEE2 then 
1^ Generate the tiles ip, t\ and put it to Tg^ 

else if tiletype — TWEl then 
1^ Generate the tiles if), t'i,t'2,t'^, t\, t'^ and put it to Tg^ 

else if tiletype = TWE2 or tiletype = TEE2 then 
1^ Generate the tiles t'Q,t\ and put it to Tg^ 

Generate the tiles in the Encoding Area and put it to Tg^ 

if tiletype = DEEl or tiletype =^ DEE2 or tiletype = TEE2 then 
1^ Generate the tiles between glue A and glue gi^, put those tiles to Tg^ ■^^\ 

else if tiletype = TWEl or TWE2 then 
1^ Generate the tile with glue gin, put it to Tg^ 

return T^^^^,; 



Algorithm 13: Generate-connection-tile() 
Input: tiletype, the type of the tile 
gin, the input glue 
gout, the output glue 
e„, the code of the north glue 

maxbits, the number of binary bits to encode all of the glues G Gns 
Output: Tg^ j^j, Tile set at temperature 1 in 3D 

nd,it^<i>-' 

if tiletype = SWN2 or tiletype = SEN2 then 
|_ The number of tiles to be generated at the bottom of dotted box in the figure is (maxbits x 2). 

else if tiletype = DWN2 or tiletype = DEN2 then 

Generate the tiles showed as the figures in the Table[3] with the input glue gin and output glue gout- 
Put all of the tiles to Tg^^^. 

return TL^ 



37 



Figure 9: The logic binary tree for constructing decoding tile set (direction left). 



Algorithm 14: Generate-encode-tile-fixed() 
Input: tiletype, the type of the tile 
gin, the input glue 
gout, the output glue 
e„, the code of the north glue 

maxbits, the number of binary bits to encode all of the glues E G. 
Output: Tg^ j^j, Tile set at temperature 1 in 3D 

Generate the tiles with the input glue gin, 

Generate the tiles in the Encoding Area and put it to Tg^ 

Generate the tiles with the output glue gout, 

return T^j , , ; 




Figure 10: The logic binary tree for constructing decoding tile set (direction right). 
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Figure 12: The tile set to decode one bit of the code (Direction right, K=4). 
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